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Abstract 

Using the Closed-Time-Path approach, we perform a systematic leading order 
calculation of the relaxation rate of flavour correlations of left-handed Standard 
Model leptons. This quantity is of pivotal relevance for flavoured Leptogenesis 
in the Early Universe, and we find it to be 5.19 x W'^T at T = 10^ GeV and 
4.83 X IQ-^T at T = lO^^GeV. These values apply to the Standard Model with a 
Higgs-boson mass of 125 GeV. The dependence of the numerical coefficient on the 
temperature T is due to the renormalisation group running. The leading linear and 
logarithmic dependencies of the flavour relaxation rate on the gauge and top-quark 
couplings are extracted, such that the results presented in this work can readily be 
applied to extensions of the Standard Model. We also derive the production rate of 
light (compared to the temperature) sterile right-handed neutrinos, a calculation 
that relies on the same methods. We confirm most details of earlier results, but 
find a substantially larger contribution from the t-channel exchange of fermions. 
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1 Introduction 



The origin of the matter anti-matter asymmetry in the Universe is one of the most 
important open questions in Particle Physics and Cosmology. Leptogenesis [T]!^ is among 
the most plausible possible mechanisms for generating the observed baryon asymmetry 
dynamically in the Early Universe. The recent discovery of a Standard Model (SM) 
like Higgs boson provides further support for this mechanism, where the asymmetry is 
generated by out-of-equilibrium decays of heavy sterile right-handed neutrinos into a 
Higgs boson and a lepton. 

In recent years, substantial progress was made in the theoretical description of Lepto- 
genesis using Non-Equilibrium Quantum Field Theory in the Closed Time Path (CTP) 
formalism [SI [6l [7]. Particular attention was given to the derivation of quantum evo- 
lution equations for distribution functions [H |9l [101 El 112], the calculation of thermal 
corrections to the charge-parity (CP) asymmetry [131 EH [151 [HI [13 [ISl [IH], and related 

questions [2ni[2Il[221[2Sl[2aESl[2Sl[2Z]. 

Interaction rates mediated by Higgs- Yukawa couplings are crucial for Leptogenesis 
predictions. The production and decay rates of the lightest right-handed neutrino de- 
termine the strength of the washout, and they are mediated by their Yukawa couplings 
Y to SM lepton doublets and Higgs bosons, while the lepton flavour equilibration rate 
determines the temperature scale where flavour effects become relevant for calculations 
of the asymmetry [281 [SSI [30] , and it is mediated by the SM lepton Yukawa couplings h. 

At low temperatures, both processes [right-handed neutrino (inverse) decays and 
flavour equilibration] are described at leading order (LO) through 1 o 2 processes, 
namely decays and inverse decays, involving only the Higgs boson and left- and right- 
handed leptons and sterile neutrinos. At high temperatures, both rates are phase space 
suppressed. Therefore a complete LO calculation of these crucial interactions must 
include processes where additional gauge bosons are emitted or absorbed or where an 
intermediate Higgs boson decays into a pair of top quarks. For right-handed neutrino 
(inverse) decays, the LO analysis pertinent to high temperatures has first been performed 
in Refs. [211 [32] • In some earlier and also more recent papers [2^1 [331 [311 [3S] , the inclusion 
of the gauge interactions is restricted to the effect of the modified thermal dispersion 
relations (thermal masses), but these approaches should only partly capture the LO 
effects. Closely related to the problem of right-handed neutrino production and the 
relaxation of the flavour correlations of active leptons is the production of photons in a 
quark-gluon plasma, which is calculated e.g. in Refs. [3S1[3Z]- 

Gauge boson corrections in the low temperature regime have recently received some 
attention in the works [381 [39] ? where a non-relativistic analysis including leading rela- 
tivistic corrections has been performed. When the right-handed neutrino mass cannot 
be neglected, the tree-level scattering amplitudes exhibit infrared (IR) divergences, and 
the familiar cancellation of these with contributions from virtual diagrams has to be 
generalised to a finite-temperature environment. It turns out that the cancellation of IR 
divergences in the non-relativistic regime can be generalised to fully relativistic processes 

^See e.g. Refs. [ll[31ll] for recent reviews. 
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as well [10]. 

The main topic of this paper is the calculation of Higgs Yukawa mediated inter- 
action rates in a high temperature background. We have performed this calculation 
using methods based on the two-particle-irreducible (2PI) CTP formulation of Non- 
Equilibrium Quantum Field Theory, within the framework that has been applied to 
calculate CP- violating rates in Refs. [HI [19]. A main advantage of the 2PI approach is 
that it directly gives rise to a derivation of the appropriate description of the physical 
screening that regulates the t-channel divergences that occur in certain tree-level scat- 
tering diagrams. The screening is implemented by the use of a resummed propagator, 
which in contrast appears as an ad hoc prescription within an approach based on the 
calculation of scattering matrix elements. 

Presently, our main phenomenological interest is to perform a LO calculation of the 
damping rate of lepton flavour coherence, that is of importance for flavoured Leptoge- 
nesis. This flavour damping rate can be used as an input to the systematic analysis of 
flavour decoherence in Leptogenesis, that has been presented in Ref. [TT]. Having ob- 
tained the various diagrammatic contributions to flavour damping, it is a simple matter 
to extract the (inverse) decay rates of right-handed neutrinos in the high temperature 
regime. We thereby confirm most details of the results presented in Refs. [311 [32], but 
find a substantially larger contribution from t-channel fermion exchange. The latter dis- 
crepancy should be due to a different technical implementation of the extraction of the 
coefficients of the logarithmically enhanced (in squares of the gauge coupling constant) 
and the standard perturbative contributions. In order to test our method, we have there- 
fore compared its result to a fit to the numerically performed integral and found it to 
agree within the expected accuracy of the present approximations. 

The remainder of this paper is organised as follows: In Section [21 we introduce 
the formalism and explain how the relevant interaction rates are obtained in this setup. 
Contributions to the interaction rates from self-energy type and vertex type diagrams are 
calculated in Sections [3] and [4] respectively, while coUinearly enhanced 1—7-2 processes are 
calculated in Section [51 The final numerical results and implications for phenomenology 
are presented in Section [SI before we conclude in Section [71 

2 Setup 

2.1 Goal of the Calculation 

The interactions that lead to the production of right-handed neutrinos Ni are the Yukawa 
couplings Yia- The lepton asymmetry can be partly transferred from left-handed lepton 
doublets ib to right handed active leptons Rq through the Yukawa couplings hab- These 
also lead to the decoherence of off-diagonal correlations between the lepton flavours. 
Explicitly, these Yukawa interactions are given by the Lagrangian terms 

-Yiai^Ni4>^Pl,i:ia - hab<P^KaPhiJlb + C.C. , (l) 
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where c.c. stands for complex conjugation. Four-component spinors are denoted by ijj, 
and the subscripts indicate the fields that they are associated with. While a formulation 
in terms of Weyl spinors is also possible and perhaps more appropriate, we choose here 
four component spinors in order to make use of the standard identities for Dirac matrices 
for the simplification of the spinor algebra. The indices of SU(2)l are suppressed, the 
contraction of these between the fields and i is implied, and = (e0)^, where e is the 
antisymmetric rank- two SU(2)l tensor. 

The goal of this paper is to calculate interaction rates mediated by Yukawa couplings 
in a finite temperature, finite density environment like the Early Universe. Tree level 
1^2 processes that are proportional to are kinematically suppressed in the Early 
Universe, since all involved particles are massless {i.e. masses much smaller than the 
temperature) at temperatures T above the Electroweak phase transition. Similarly, when 
T ^ Mat, where is the mass of the right-handed neutrino, the right-handed neutrinos 
can be approximated as massless, such that decays N (pi are suppressed. 

With the tree-level 1 <(-)■ 2 channels forbidden or suppressed, the LO rates involve the 
radiation of one additional gauge boson, which leads to scattering rates that are para- 
metrically of order g'^h'^ or g'^Y'^, respectively, where g denotes a gauge coupling. At the 
same order, there are also contributions from collinearly enhanced 1^2 processes when 
medium effects mediated by gauge interactions are resummed into the propagators |31j . 
Due to the sizable top-quark Yukawa coupling, a leading order calculation should also 
account for the Higgs boson decaying into a top-quark pair. 

For the present calculation, we assume that all external particles are massless (be- 
fore including thermal corrections), what leads to significant simplifications {i.e. the 
vanishing of numerous contributions from the CTP Feynman rules and the absence of 
soft and collinear IR divergences) that become apparent during the course of the calcu- 
lation. Within the CTP formalism, the production and relaxation rates can be inferred 
from the collision term, as explained in Section 12.31 The collision terms for A^, i and 
R encompass similar diagrams, that are related among one another by an exchange of 
the Yukawa coupling matrices and the gauge coupling constants. In the following, we 
therefore derive the collision term for i, but we factorise the results in such a way that 
they can easily be employed in order to obtain the production and destruction rates for 
R and as well. 

2.2 Definitions 

We perform the calculation of the collision terms within the framework developed in 
Refs. [m [19]. Here we just give the definitions that are relevant for the present calcula- 
tion. The fermionic Wightman functions that appear in the collision term are 

iS<{k) = -2S-^{k)f{k) , (2a) 

iS>{k) = 2S^{k)[l-f{k)], (2b) 

where S'^{k) is the spectral function that determines the location of the quasi-particle 
poles and f{k) is the distribution function that will be specified more precisely below. 
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Tree- level propagators, which we indicate by a superscript (0), are recovered when re- 
placing the spectral function S'^ with a delta distribution: 

S^^^^k) = ttPx (I^ + m)6{e - m2)sign(A;°) . (3) 

In case of massless, chiral fermions, either Px = Pl,r (^-C- Pl for i and Pa for R); 
whereas for Dirac and four-component Majorana fermions, Px = 1- The relations for 
the distribution functions are 

„,, f/(k) for i»>0 „ ■ . 

^(*) = |l-/(k) to. *«<0 

where /(k) and /(k) are the distribution functions of particles and anti-particles. In 
kinetic equilibrium, the distributions are of the Fermi-Dirac form 

where /3 = 1/T. Note that and f{k) may be hermitian matrices in order to take 
account of flavour coherence [TT] . 

To calculate the self energies, we also need the bosonic Wightman functions for the 
Higgs and gauge bosons. For scalars, these are 

iA<(A;) = 2A-^(A;)/(A;), (6a) 

iA>(A;) = 2A-^(A;)[1 + /(A;)], (6b) 

where the tree-level spectral function for scalar particles of mass m is given by 

A(°)'4(A;) = 7r5(A;2 - m')sign(fc°) . (7) 

The distribution function f{k) with the four-momentum as an argument can be related 
to the particle and anti-particle distributions /(k) and /(k) as 

,,,, f/(k) for A;0>0 

/W = jlV+/(k)) fo. .»<o 

For distributions in kinetic equilibrium, 

^(^) = - 1 • (9) 

Finally we use gauge-boson propagators in Feynman gauge, such that 

iA<^> = -(7,.iA<'> . 

Evolution equations for the distribution functions f{k) can be derived from the 
Dyson- Schwinger equations for the Wightman functions S'^^ and A"^. In Wigner space 
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Figure 1: Diagram that represents the flavour- sensitive contribution to the lepton self 
energy that occurs at one-loop order in the 2PI approach. Double lines indicate 
exact (in the present approach one-loop resummed) propagators and the external legs 
are understood to be amputated. 



and to flrst order in gradients, the relevant equations for the the lepton doublet and for 
the right-handed neutrino are 

idn'Sf> - [k . 77° + ?f 7°, ^l'Sr>] - [z?<'>/, I'Sf] = -\ {^C, + zCj) , (10) 

ida'^s^^^ = -iCn. (11) 

For a derivation of these equations and the further treatment of the left-hand sides we 
refer the reader to Refs. [HI [19]. Physical processes, like the Yukawa induced scattering 
and decay processes that we are interested in here, are encoded in the collision terms C, 
as explained in more detail in the next Section. 



2.3 Contributions to the Collision Term 

In this Section and throughout most of the remainder of this paper, we consider the 
flavour relaxation rate for left-handed SM leptons I. The methods employed and the 
calculations performed readily give rise to the production rate of sterile right-handed 
neutrinos A^, as it is explained in Section 12.51 Numerical results of this production 
rate together with a comparison to the earlier results of Refs. [311 [32] are presented in 
Section 16.21 

The collision term for the lepton doublets, 

C,{p) = i?>(p)i5<(p) - i^f{p)iS>{p) , (12) 

can be calculated perturbatively by employing a loop expansion of the leptonic self en- 
ergy i^f'^- Furthermore, the propagator iS*/'^ contains self-energy insertions that may 
be of importance and that can also be expanded perturbatively. In this paper, we restrict 
ourselves to flavour-sensitive contributions that are at least of order h"^ in the Yukawa 
couplings and that we indicate by a superscript fl, i.e. and Cf. Due to the non-trivial 
flavour structure, these couplings will give rise to flavour decoherence, in contrast to 
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the flavour- diagonal gauge couplings. The latter interactions are however of importance, 
as they open up the phase space for flavour-decohering processes (in combination with 
the coupling h) that would be kinematically forbidden otherwise. Moreover, gauge in- 
teractions maintain kinetic equilibrium for gauged particles, which is what we assume 
throughout this work. 

Up to two loops, there are contributions to the flavour-sensitive rates from the one- 
loop self energy in Figure [1] and from the two-loop vertex-type diagrams in Figure |5l 
which are obtained from the self-energy diagram by connecting two different propagators 
with a gauge boson. 

Note that in the 2PI approach, no two-particle-reducible diagrams appear where a 
gauge boson connects to both ends of the same propagator. Similarly, there is no dia- 
grammatic contribution from the insertion of a top-loop into the Higgs boson propagator. 
While such diagrams with a self-energy insertion do not derive explicitly from the 2PI 
loop expansion of the effective action, self-energy insertions are of leading importance 
and are accounted for implicitly, by using resummed propagators to evaluate the collision 
terms. In fact, self-energy insertions in all propagators present in Figure [T] contribute to 
the flavour decoherence at LO. As we will see below, it is possible to expand the Higgs 
propagator up to the order of single loop insertions, whereas the lepton propagators must 
be maintained in the fully resummed form due to the presence of t-channel divergences 
from fermion exchange. 

According to the loop expansion indicated in Figures [T] and El we decompose the 
collision term as 

cf = cf^ + c^'^*"^ (13) 

The contribution from the lepton self-energy diagram. Figure [1], is given by 

itf^^'ip) =h^ J 0^,^A^^mP - k - q)iAf{k)iS^\q)h , (14) 

i.e. 

= ifP'-{p)iS<{p) - ifp<{p)\S>{p) . (15) 

Within the 2PI approach, the propagators iS'^,R and iA,^ are understood as the exact 
propagators. In the present work, we approximate these by propagators that contain the 
resummed one-loop corrections that arise from gauge interactions and from top-quark 
loops, cf. Section [3l Diagrammatically, the coUison term f|T5|) is given by Figure |2l where 
the double-hne propagators represent the resummed propagators. 

When all propagators are on shell and massless, = /c^ = = 0, this self energy 
does not contribute to the collision term, so that no contribution of order h'^h arises. 
However once the medium effects are included by replacing the tree level propagators 
with resummed propagators, one obtains contributions that are parametrically of order 
g'^h^h. The LO terms of this type will be calculated in Section [31 
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Figure 2: The collision term Cf^\ Eq. (pjl. 

The vertex-type diagrams give an order g^h)h contribution to the collision term when 
evaluated using tree-level propagators. While in principle, it would be possible to eval- 
uate also these diagrams with finite temperature resummed propagators, the difference 
between the two calculations is of higher order in the gauge coupling. For the purpose of 
the present work, it is therefore sufficient to calculate C™'"*'^'^ using tree level propagators, 
which will be performed in Section HI 

In Ref. [31] it is shown that coUinearly enhanced \ ■k^ 2 processes appear as a 
third contribution that is parametrically as important as the two contributions discussed 
above, and which requires a resummation of ladder diagrams with an arbitrary number 
of soft gauge bosons inserted between the Higgs and the lepton propagator in Figure [1] 
These multiple scattering contributions will be evaluated in Section [51 



2.4 Flavour Equilibration Rate 

The evolution equations for the left- (right-) handed lepton asymmetries qi {q^) are given 

by [m 

{W,, qt\ + 2S, - Y\ , (16a) 

(16b) 

where in general, and Wi^ S^ and F^j^ are 3x3 matrices in fiavour space. The 
charge densities g^^R are matrix-valued generalisations of charge densities of left- and 
right-handed SM leptons, that can thereby also account for fiavour coherence. The 
lepton-number violating washout term is denoted by Wg^ the source term for the CP- 
asymmetries by S^, and F^j^ describe the equilibration of lepton number between I and 
R as well as the decay of fiavour coherence. In the standard scenarios of Leptogenesis, 
there is no direct source for a charge asymmetry in right-handed leptons R and no direct 
lepton-number violating coupling that woul induce a washout. 

Here we are interested in the fiavour equilibration rates F^j^ that depend on the 
Yukawa couplings hah- The rate F^ can be obtained by integrating the relevant part of 



It 
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the collision term Cf. 

74^ 1 r 



+ h.c. 

(17) 



where h.c. denotes hermitian conjugation in flavour space, and the superscript fl indicates 
that we only consider those contributions to the collision term that are at least second 
order in the charged lepton Yukawas hab- 

We assume that the particles £ and R are in kinetic equilibrium, i.e. their distri- 
butions are of the Fermi-Dirac form with generalised, matrix-valued chemical potentials 
f^£,Rab- Moreover, we assume that the chemical potentials are small, fii^Rab T. We can 
then expand the collision term to linear order in the chemical potentials using 

feab{k) = ^pj^ ^ Sab + Sfeab{k) , (18) 

^f^-b{k) = /3pfcT^Yp/i£a6 = QP'j^^k^r^Qiab • (19) 

Note that this also implies a decomposition of the Wigner functions iS^'^ = iS^J,^ + i5Si, 
where the deviation from equilibrium is given by 

i6Se = -2Sf6fe{p) . (20) 

Accordingly, the term linear in the right-handed SM lepton asymmetry is obtained 
by expanding iS^'^ , which enters through Eqs. (IMI) and (ITTI) . Now we expand 
Eq. ( fTTI) in fii^^i/T. Due to the Kubo-Martin-Schwinger (KMS) relations, iS^^^^{p) = 

—G^^°iSf^^^{p) and i^^^^{p) = — e^^°i|^^^^q(p), where i?^^^eq is the contribution from the 
equilibrium parts of the individual propagators to the self energy, the leading terms in 
this expansion are linear in [i.e. the equilibrium contributions vanish). Therefore, 
the flavour- sensitive contribution to the lepton collision term takes the general form 

J use MR 

= V (^^^^^ + ^•"•) + V ^^^^^^ + ^-'-^ ' ^^^^ 

where we have defined the reduced scattering rates ■y^^^'^. Due to lepton number con- 
servation of the interactions mediated by h in combination with flavour-blind gauge 
interactions, we expect that 7^''^ = —'-f^^^^ what we explicitly verify in the calcula- 
tions. Note that qiab = ^^tab ~ ^^7ab is difference of the deviations of the lepton and 
anti-lepton number densities from thermal equilibrium. Fast gauge interactions ensure 
that fn] 5n^^^ = —SnJ^fj such that F^j^ can be obtained from the collision term according 
to Eq. (fT7j) . otherwise the collision terms for particles and anti-particles would have to 
be evaluated separately, and equations of motion for each ^n^^ and would have to 
be solved. 
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The coefficients ■y^^^'^^ have contributions from several diagrams, 

^RSiM ^ ^fl(A)5£,5R _^ ^fi(B)5£,5R _^ _ _ _ ^ (22) 

which we calculate in the following. Likewise, we decompose the self energies the 
collison terms Cf and the damping rates into various diagrammatic components. 

2.5 Right-Handed Neutrino Production Rate 

The kinetic equation for the right-handed neutrino density JnAp) is given by 

|/^.(P) = A(P), (23) 

where the decay (or inverse decay) rate Di is obtained by integrating the collision term 
over p° [IS] : 

1 fdp\:^.^r^O, 



A(p) = -^signipyrCmip) (24) 

= 11 ^sign(p°)tr \irNuip)^SUp) -^rm^iP)^SUp)] ■ 

Here, we neglect the effect of possible off-diagonal correlations in the neutrino propagator. 
As long as the mass differences M^rj — Mjsfj are much larger than the relaxation rates D, 
this is a suitable approximation, because the off-diagonal correlations oscillate rapidly 
and do not give a coherent contribution to the diagonal evolution. In the case of a strong 
mass-degeneracy, the full evolution of diagonal distributions as well as the oscillatory off- 
diagonal correlations must be considered, as it is described in Ref. [21]. For simplicity, 
we suppress in the following the flavour indices for the sterile right-handed neutrinos. 

For the right-handed neutrinos, the medium corrections to the propagator are pro- 
portional to \Y\^ and therefore negligible in general, such that we can use the tree prop- 
agators ([3]) with mass in the Wightman functions ([2]). Inserting these into Eq. ( !23|) . 
the integral can be performed analytically, and one obtains (suppressing the flavour 
indices) 



^^Nip) (::^^^(i-Mp))-(i-^m?)Tt'^^^p^ 



(25) 

where p° = uj{p) = ^^/^~l^T4j^ and where we have used that i and are in thermal 



equilibrium, such that the self energies satisfy the KMS condition i^^(p) = — e^^°i^^(p). 



Furthermore, the spectral part of the self energy is defined as 

1 



2 (i?^ - i?^) • (26) 
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Taking account of initial conditions for fj^n and of the expansion of the Universe, solving 
Eq. fl25|) results in a non-equilibrium distribution function /Ar(p), which can be computed 
numerically, but for which no simple analytical form applies in general. Presenting an 
illustrative result for the neutrino production rate therefore is to some extent a matter 
of definition. For the ease of comparison, we follow Ref . [21] , and choose the production 
rate for /Ar(p) = 0. Approximately, this is the relevant rate for the production of singlet 
neutrinos in the weak washout scenario, provided it is assumed that initially, the density 
of these particles vanishes. The differential production rate for the sum of the two spin 
orientations is 

where 

i^< (p) = -2-^J^ti{p) . (28) 

To see more explicitly how this rate is related to the flavour relaxation rate defined in the 
previous section, it is instructive to consider the leptonic collision term Cf{p). As before, 
we expand 15"/'^ = 15"/^^ + iSSi and recall that i6Si can be approximated as being linear 
in q^. To leading order in deviations from equilibrium, we can therefore assume that the 
self energies satisfy the KMS relation. This part of the collision term then simplifies to 



1 1 = -J^i^'iP^^tApWiip') + 0{6fn) + h.c. , (29) 



with p = ±|p|. This has the same structure as the equation for the differential right- 



handed neutrino production rate ( l25l) except for the different statistical weight functions. 
A relation between tr[|i|^^] and tr[|i|!l^'^] is obtained by isolating the dependence on the 



coupling constants. In most cases, these are simple proportionalities, whereas for t- 
channel fermion exchange, there is an additional logarithmic dependence on the gauge 
couplings which can be isolated as well. Therefore, once the differential (in |p|) flavour 
relaxation rate is known, the production and decay rates can be obtained by simply 
performing the integral with the appropriate statistical weight and a rescaling of the 
couplings. 

3 Self-Energy Type Contributions 
3.1 One-Loop Self Energies 

The resummed form of the spectral function of a massless chiral fermion is [TTl UH |32] 

2 (^1^ - ^) S-^ ■ (A; - S^) - ^ (^^ - + 



S-^ik) = Px- 



4 ■ (A: - S^)]' 



(30) 
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and of a massless scalar 



(31) 



We express the spin-^ fermionic self energies as ^ = In the approximation of 

massless particles in the loop, the spectral self energy for leptons £ or R is given by 

GT2 



{k) 



167r|k| 

GT^ 
IGyrlkl 











M) 








|k| 









0\2 



2|k|T 



-/o 



(k^ 






n 



where 

/o(y°,y) = -%^-(y°)V-l/ + log 



ea 



l + e^to"-^) 



log 



1 _ qW+v) 



vr 



(32a) 
(32b) 

(33a) 
(33b) 



+Re 



a;(log(l + e") - log(l - e^"^ )) + Li2(-e") - Li2(e 



x—y" 



For the case of the Standard Model lepton doublet £, G = ^{3gl + g'f), whereas for 



right-handed leptons R, G = 2gl. Notice that these values for G include a factor of 
two that accounts for the polarisation states of the gauge bosons, as the spectral self 
energies are defined here for a single bosonic and a fermionic (with both spin states) 
degree of freedom in the loop. Note as well that we distinguish between the lepton self 
energy (l32l) . that is flavour-diagonal and of order g'^, and the flavour- sensitive self energy 
itfi, for which the LO terms are ~ h^g^ and ~ h^g^logg^. 
The spectral self energy for the Higgs field is given by 



_ 2^2 + 2^1 


,r, 




2, 

^log- 


IGtt 








ki - 


^log 


1 + e^ 


167r|k| I 


l + e'^ 



2 



e*^ 2 



(34a) 



for >0, 



|k| + fc" 



IGtt |k| \ /3 
3/1? e / 2 1^ 



1 -e' 

3 |k| + fc° 



|k|-fcO 



(34b) 



for A;^ < . 
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Figure 3: Diagrammatic representation of the LO expansion of the one-loop 2PI coUision 
term, Eq. flgT]) . 

When we substitute tree-level propagators for all three fields R and that appear 
in the expression for C|°'^ defined by Eqs. ( I14|[T^ . we obtain a vanishing result, because 
the 1^2 process is kinematically forbidden (when neglecting the tree-level masses) 
without including the finite-temperature corrections. At leading order in the gauge and 
top-quark Yukawa couplings, the medium corrections add linearly such that the one-loop 
(in the 2PI sense) self energy f[T^ for the lepton doublet can be approximated as 



i^f(p) 



for = 



where 



(27r)4 (27r)4 

d'^k d'^q 
(27r)4 (27r)4 

d'^k d^q 



{2nY6\p-k-q)i^f''\k)\S^^{q)h, 
{27,r5\p-k-q)iAf{k)iS^;^-\q)h, 
{2^r6\p-k-q)iAf''\k)iS^^^^''\q)h. 



(35) 

(36a) 
(36b) 
(36c) 



(27r)4(27r)4' 

Here, we have also defined the contribution from on-shell Higgs bosons and right-handed 
SM leptons i^^°'^^'^^(p), where it^f^^^'^{p) = for = for the kinematic reasons 
mentioned above. However, since the resummed propagator iS^'^{p) is non- vanishing 

for 7^ 0, there occurs a contribution involving i'^^^^^"''' (p) that is kinematically allowed 
due to gauge bosons that may radiate from i. In summary, we can decompose the 
self-energy contributions to the flavour- decohering collision term as 

^ Cf^^ + + C?(^) , (37) 



where 



Cf^-^^ =tr 



i#'^^>(p)i5f <(p) - i4^'^)<^^^^c(o)>^ 



{p)iSr'{p) 



C 



=tr 



\t'r'^{p)is<{p)-itr{p)is>{p) 



(os)<. 



(38) 
(39) 



A diagrammatic representation of the decomposition (IS7I) is given by Figure O 
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3.2 Scatterings via the Higgs Boson 

Substituting the tree-level spectral function ([3]) for I and R and the one-loop resummed 
spectral function for 0, the collision term becomes 



^ (27r)4(27r)4^ ^ A;^ + [^^(A;)]2■ 

x 27r(5 (p^ ) 27r 5 ( ) sign (p° ) sign ( g° ) 

x/^t [(1 + f^ik'))Mq')hil - - /<^(A:°)(1 - /R(g°))/i/,(/)] . 

It is now useful to notice that n-^(/c) is first order in k"^ for /c^ — [c/. Eq. (134|) ]. as well 
as 2p ■ q = k"^ (due to the relations imposed by the (5-functions). To leading order in the 
couplings and hf, we may therefore replace 

2nrf 2nf ^^^^ 



*-4 + |n^]2 fci • 

such that in this approximation, C^^'^^ is proportional to a-nd ^f- Alternatively, 
one can derive Eq. (l40l) with the replacement (14T!) from a CTP two-loop (two particle 
reducible) self energy in terms of tree-level propagators. 

We now assume that 0, i, R are in kinetic equilibrium, and that the chemical poten- 
tials of i and R are small compared to the temperature, such that we may relate 

where 5fi^^{jP) = fi.B^ijP) — fp^iv^)- Note that 5fi^^{jP) and are understood as ma- 
trices in flavour space and fp^ip^), the Fermi-Dirac equilibrium distribution for vanishing 
chemical potential, is therefore implied to be proportional to the unit matrix. 
With the definitions ( IT71l21ll22p . it follows 

d'k f d^p f A.^2vr)V(p + A:-g)^^^^27r5(p2)27r(5(g2) (43) 



(27r)4 J (27r)4 J (27r) 



X sign (p° ) sign (g° ) — p „ 6/3 



=7.71 X 10-^ X {^gl + ]^gi^ T + 1.32 x IQ-^ x h^T , 
where the last expression is the result of numerical evaluation of the integrals. 
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3.3 Scatterings via Fermions 

Now consider the terms that arise from substituting the tree-level spectral function ([7]) 
for the scalar propagator and ([3]) for one of the fermions ^ or R, while using the resummed 
spectral function ( 130|) for the other fermion. For definiteness, we calculate the coUison 
term for scatterings via R first [i.e. with iS^ as in Eq. ( 130|) ]: 



(27r)V(p- A; + g)2tr [S^ 



(44) 



(27r)4 (27r)4 ' 
X 27r5(j9^)27r5(g^)sign(p°)sign(g°) 

X [h^h{k')h{l - /,(/))(! + /^(g°)) - h\l - fK{q'))hMp')U{q')] . 

The contribution from scatterings via i can directly be inferred from the evaluation of 
this term. Expanding in the deviations of from chemical equilibrium, we write 

J-(fc°,p°,g°) = [h^Mk')h{l - /,(/))(! + U{q')) - h\l - h{q'))hf,{p')U{q')\ (45) 

1 1 



/itM/,(/) 



+ 



+h^5Uk'')h 1- 



+ 



1 



e/3p« + 1 Q^q" _ 1 



what leads to the decomposition C^^^"* = Q^'^^^^^ _|_^fl(R)5R accordingly for and 7^. 

The collision term (jUj) can straightforwardly be evaluated numerically. After making 
use of the 5-functions, homogeneity and isotropy, two numerical integrations are left. 
However, it is useful and instructive to isolate the dependence on the coupling G. 

It is known that the phase space integrals over the scattering matrix elements exhibit 
a logarithmic divergence for zero-momentum exchange of a lepton R in the t channel. 
The scattering approximation is recovered in the present approach when omitting the 
self energies in the denominator of the spectral function, Eq. (130|) . which would lead to a 
logarithmic divergence in the integral (1441) for k^ = |k| =0 and /c^ < 0. A simplification 
corresponding to the replacement 0411) is therefore not suitable for the present integral. 
A finite part can however be extracted when subtracting those terms from the spectral 
self energy that are not vanishing in the limit /c*^, |k| — )■ 0. These are precisely the 
contributions that are of the form of the hard thermal loop (HTL) approximation and 
that we indicate by a tilde. We therefore define 



vr 



(46) 
(47) 



The barred quantities S'^ and E are defined through the replacement Ii — Ji, and 
accordingly S-^ = S-^ + t^, cf'^'^ = cf"^^ +C'~' 
approximation. 



fi(R) 



Note that Jq vanishes in the HTL 
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For the purpose of calculating C^^^^ to leading order in G, it is sufficient to approxi- 
mate 



(48) 



The result for Cf^^^ is therefore manifestly proportional to G. 

It remains to calculate the part that originates from the term of the HTL form. 
Therefore, we need to find 



=27r(5(p2)sign(p°) 



(27r)4 (27r)4 



(27r)^(5(p -k + g)27r5(g^)sign(g°) 



(49) 



x2tr 



m - ■ (k - ) - i-" 



R 



2 ^^3 



i ^R 



4 



In order to extract the dependence on G, it is useful to split the integral in a region 
where |k| > k^ ^ VGT, where the denominator simplifies (because the self energies may 
be neglected there far from the single-particle poles) and a region where |k| < fc* <C |p|, 
where the angular integration simplifies, C^^^\p) = C^^^^^{p) + ^^^>fc. (p)- Notice that as 
a consequence of this split, the result should only be valid when |p| ^ y/ GT. Due to the 
phase space suppression, this does however not spoil the LO calculation of the flavour 
relaxation rate. 

Dropping the terms oc G in the denominator, we evaluate 



^fl(R) 



GT^ 

fi=\p\ 2'^7T 



dcos'd 



d\k\ k^ 


- 2 


Pl 


lkl 




p 


— 





(50) 



where 



IpI — a/p^ + — 2|p| |k| cos'd . 



(51) 



(Note that the HTL contributions are only present for k"^ < 0.) The collision term 
for < may be obtained when noting that it is even in provided the particle 
distributions are in kinetic equilibrium. 

The logarithmic dependence can be isolated through integration by parts. We sim- 
plify the angular integration and additional terms using /c* ^ |p|, such that k^ ~ 
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k| cos'd, and obtain 

fl(R) 
k 



pnyrx) , ^fl(R) 

'-'^,>fe,LOG "I" W,>A:*FIN ' 



c 



fl{R) 
,>fc,LOG 



^fl{R) 
W,>fc.FIN 



2% 



27r5(p^)log(/3A;,)^(0,p^ V 



(52a) 
(52b) 



27r(5(p2 



dcos'& j rf|k|log(/3|k|) 

^fl(R) 



d 



2|p| 



a(/3|k|) fco- IpI 



(52c) 



which we have separated into a contribution Cg >fc^LOG ^'^^^ depends logarithmically on 

and an integral C^l^^^pjj^ that depends only linearly on [i.e. that is finite for k^ — ^ 0), 
such that we can take the lower bound of the integration to zero when k^, <^ |p|. 

In order to calculate X<p^ , it is necessary to take account of the screening that is in- 
duced by the self energies. In addition to the spectral self energy, also the hermitian part 
is of importance. Because -C T, it is sufficient to consider the HTL approximations 



j:"\k) 



GT^ 



log 



32|k| 

GT^k^k' 
32|k|3 





|k| 


A;0- 


|k| 



log 





|k| 


GT^k' 


A;0 - 


|k| 


16k2 



(53a) 
(53b) 



We approximate |p — k| = |p| — p ■ k, and moreover, we evaluate the statistical 
functions for k'^ = |k| ^ 0, as it is appropriate for k^ <^ |p|. Numerically, we can then 
obtain the value of 



r'm) 



' (27r)2|p| 



dcos^ j k"c/|k|tr pS^{k) 

1 



(54) 



The prime on S^{k) indicates that we evaluate this expression by replacing G — ?■ G', 
and the tilde indicates, that we use the HTL approximation. Above expression then 
corresponds to the infrared contribution to the scattering rates that arises from a UV 
cutoff k^ and a squared coupling G'. From the dependence of the HTL self energies on 
G and on the four-momentum, we observe that a simultaneous rescaling of the squared 
coupling by a factor of G/G' and of the momentum by ^fCjC' rescales the value of 
the integrand (including the integration measure) by an overall factor of G/G' . There- 
fore, [G / G')C'^^'^ also describes scatterings for a squared coupling G and a UV cutoff 
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•\/ GjG'k^. In order to obtain the contribution for a cutoff A;*, we must add the integral 



^TfT^,, , , = -'^^sip T r ' i ^cos^ / k^^ikitr ps^ik) (55) 



^^(0,pO,-pO 



^e[y/G/G^K,k,] ^ (27r)2|p| 



In summary, the wave-function contribution to the scattering rate can be decomposed 

as 

^fi(R) _ ^fi(R) ^fl(R) ^fi(R) ^fl(R) _^^,fl(R) , , 

The terms C^''^'', C^^yl^^m ^^"^ ^^'i<jt' should be evaluated numerically and are pro- 
portional to G. The logarithmic dependence on G, that results from the screening of 
scattering processes with small momentum exchange is isolated in 

CT'iLoc +C«<^,,,, = 2.^(p^)|^.og {^ffk?) V) . (57) 

We emphasise that the final result Cf^^^ is by construction independent on the choice of 
k^, k'^ and G', as the dependence of above expression on these parameters is compensated 
by • Recall that this approximate cancellation of the dependence on /c^,, k'^ and G' 



is a consequence of the approximate behaviour of the integrand when GyT <^ |k| <^ |p|. 
Below, we verify this numerically in order to test the accuracy of the approximations. 

Scatterings may as well proceed via the exchange of a doublet lepton i. This con- 
tribution to the collision term is C^^^\p), Eq. (155]) . With the integration over d'^p, the 
relevant integrals are identical to those for the exchange of R. From the result for the 
scattering via R, we can therefore directly infer the contribution from exchanges of i. 

Substituting the collision term C^^'^^ into the expression for the relaxation rate ( l2Tl) 



and using the definition ( ITTl) . we find 

^mm = _^m,im = 4.40 x lO"^ xGT- 9.33 x io~^ x GTiogG , (58) 

where G = 2gf for R exchange (superscript (R)) and G = + \g'l for £ exchange (su- 
perscript (£)). Notice that there is an analytical expression for the numerical coefficient 
of the contribution that is logarithmic in the coupling constant, 

/ ^27r5(p2)sign(/)^ (- + ^ , = ^ 9.33 x lO'^T. 

(59) 

The independence of the result (39) on G and k^ is valid up to order GlogG. The 
next-to leading expressions are of order G^logG, such that for values of G ~ 0.3, one 
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Figure 4: Cf^^\p)py {6qe)27r6{p'^) for 6qn = over p° = |p| {p^ = 0) for G = 0.3. Blue: 
semi-analytical, Eq ( l56l) . Red: full integral. Blue, dashed: semi-analytical without the 
approximation (HHi) . 

should expect to yield an accuracy of about 20%, which is obviously less than what is 
suggested by the number of digits given in the numerical coefficients. Note that this 
estimate for the accuracy is very crude as it does not account for loop and phase-space 
factors. An estimate of the next-to-leading order (NLO) contribution is non-trivial, and 
a calculation of the NLO production rate of photons from the quark-gluon plasma has 
recently been reported in Ref. [13]. 

We can also extract the coefficients of the contributions that are linear and logarith- 
mic in the couplings by directly performing the integral ( l44l) for different values of the 
coupling. By this numerical fitting procedure, we find 

^mm = 3,72 X 10~^ X GT - 8.31 X 10"^ x GT log G . (60) 

This decomposition is valid for a range of G between G = 0.01 and G = 0.6 and can 
therefore also be used for the calculation of related processes in other Baryogenesis 
scenarios. The numerical difference between the results ( 1581) and ( 160|) is due to a partial 
inclusion of higher order effects that is implied when ( 144|) is integrated directly. 

In Figure HI we show a comparison of the result from the numerical integration for 
Cf^^\p) and the semi-analytic result Eq. (ISBj) . There is a very good agreement for 
large p^/T, while we observe the anticipated breakdown of the approximations when 
p° = IpI ^ V GT is not valid. 

We finally note that in the case of right-handed neutrino production, there is no 
t-channel divergence from the exchange of a left-handed SM-lepton i when 7^ 0. 
Rather, this contribution has a logarithmic divergence in M^, as shown in Ref. [ID] , 
indicating that the resummed propagator for i should be used as well when 7^ but 
Mn < T. 
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4 Vertex Type Contributions 




Figure 5: Two-loop self energies that contribute to the flavour relaxation rate, where V 
stands for both SU(2) and U(l) gauge bosons. Only (A) and (B) are allowed for SU(2) 
gauge bosons, but all five diagrams appear with U(l) gauge bosons. In the present 
approximation, we use tree-level propagators, as indicated by the use of single lines, and 
the external legs are understood to be amputated. 

The two-loop self energies that involve two charged lepton Yukawa couplings and 
that descend from 2PI vacuum graphs are shown in Figure [5l All diagrams are obtained 
from the one-loop self energy by connecting two different propagators with a gauge boson 
propagator. Note that diagrams (C), (D), and (E) only exist for the weak hypercharge 
gauge boson and therefore are of order gfh"^. 

Since we consider massless particles, virtual corrections to the vertices do not alter 
the fact that 1 O 2 processes are kinematically forbidden. Therefore, we can restrict the 
discussion to those configurations that contribute to the 2 -H- 2 scattering rates. Each 
diagram has two such contributions that are indicated by the cuts in Figure El 

In order to present our calculational method, we discuss in the following the contribu- 
tions that arise from diagram (C). The approach to calculating the remaining diagrams 
is very similar and therefore presented more briefly. As we do not need to use resummed 
propagators for the calculation of the LO contributions from the vertex diagrams, all 
propagators that are explicitly employed in this Section are understood to be tree-level. 
For notational simplicity, we omit the superscript (0) that would otherwise occur in 
a large number of instances. The contribution to the OTP self energy represented in 
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Figure 6: The cuts in the vertex diagrams (indicated by red double hues) that cor- 
respond to scatterings. Note that the remaining possible cuts correspond to virtual 
corrections to 1 f-)- 2 processes, that are kinematically forbidden for the (approximately) 
massless particles considered here. 

diagram (C) is 



t. 



,(C)a6 
'I 




V 



xiSi{p-q)i^f{q)\^f{p-k)h, 



(61) 



where Yj 
consider 



— 1 and Iw 



|. The collision term only depends 



on i^^'^. We therefore 





ky^^au (k + q-p) 



iS^'^{k)riS^^{p 
iS^-{k)YiS^HP 



q)iA~+{q)[k-p-qYiA--{p 
q)iA++{q)[k-p-qYiA^+{p 



k)iA+~{k + q-p) 
k)iA-^{k + q-p) 



+\s^-^{k)^^\si-^{p 




xh. 



(62) 
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It is useful to shift the momenta in the first and the third term, such that 
^Pe yp) =YRY^g^h^ ' ' 



(27r)4 J (27r)4 

- g)7^i^R+(A;)iA^+(p - k)[k - p - qriA^^{q)iA~~{k + q~p) 

-iS^+{k)YiS++{p - g)iA^+(g)[fc -p- qriA--{p - fc)iA+-(fc + q-p) 
-iS^-{p - q)riS^+{k)iA++{p -k)[k-p- g]'^iA-+(g)iA+-(A; + q - p) 

+\S^-^{k)Y-^St^{p - g)iA++(g)[A; -p- g]^iA^+(p - fc)iA++(fc + g - p)| 

xh. (63) 

In the collision term, the Dirac structure of i'^f'^{p) is dotted into the lepton propagator 
iSf'^{p), which provides a factor Taking the trace over Dirac indices and using that 
tr[7'^7''7''7°'] = tr[7°"7''7'^7'^], as well as the cyclicity of the trace, we arrive at the form 

itr[ptf'^ip)] = YRY,glh^J ^ J ^tr[^ 

X !^iS^~{p - g)7^i5j,+(A:)iA-+(p - k)[k - p - qYiA'^- iq)iA^r{k + q - p) 

-iS++{p - q)riS^+{k)iA-^-{p -k)[k-p- qr\A-^^{q)\A^^t-{k + q - p) 
-iS^~{p - g)7'^i^R+(fc)iA++(p -k)[k-p- qr\A-^^{q)\A^^^-{k + q-p) 

+iS+^{p - q)r^S^^{k)iA~^^{p -k)[k-p- g]'^iA++(g)iA;^++(A; + q - p) 

xh. (64) 

First note that it is not possible to put all five propagators simultaneously on shell. 
Now suppose that p — g is off shell, whereas the remaining momenta are on shell. Then, 
the second and third term cancel, as well as the first and the fourth. Note that this 
corresponds to the kinematic regime of the CP-violating source of the corresponding 
diagram for Leptogenesis [TT], where the cancellation does not occur in general, because 
of the complex conjugation of coupling constants. 

When three momenta are on shell, the first and the fourth term give a virtual cor- 
rection to the kinematically suppressed 1^2 process. The 2^2 scatterings are 
encapsulated within the second and the third term, where the two off-shell momenta are 
p — q and k — p. The numerator algebra of the self energy (l64l) is given by 

tr[^(^ -4){t-f- 0] = 8p-kp-q~ Ap'k ■ {p + 2g) . (65) 
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Here, we have used the on-shell conditions k"^ = = {p — k — q)"^ = 0, that apply 
for the second and the third term, while we yet allow for ^ 0, such that this result 
may be used for the application of calculating the production rate of massive right- 
handed neutrinos [10]. In what follows however, we set p"^ = 0. This leads to the great 
simplification that the zeros in the numerator and denominator cancel: 

^tr[^(^-^)(^-^-^)^] = {p-q)\k-p)\ (66) 

A consequence of the fact that this cancellation does not occur for 7^ is the presence 
of soft and collinear divergences in the particular real and virtual contributions to the 
vertex-type self energy for the production of massive neutrinos. It is shown in Refs. [Ml 
[3ni Sn] that upon summation of all contributions, these soft and collinear divergences 
cancel. 

The same discussion can be repeated for t^f'''^ by simply exchanging the first and 
the second CTP index on all propagators. For the integrated collision term, we then 
obtain 

/ (S^'"' ^ - I&I&I (67a) 

x27iS{{k + q - p)^)sign(p°) sign(A;°) sign(g°) sign(A;° + g° - p°) 
xe?(g°,fc° + g°-p°,A:°,/), 
g{E^, E2, Es, E,) = (1 + U{E,))fA{E2)h\l - U{Es))hME,) (67b) 
- /<^(i?i)(l + fAiE^m^MEsMi - ME,)) . 
Using the definition (IT7|) . we obtain the flavour-sensitive rate 

rf^) =Y^Y,gll ^/^/ ^,2^S{p')2nS{e)2nS{q')2n6{{k + q-p)') 
xsign(p°) sign(/c°) sign(g°) sign(/c° + g° — p^) 

X i (^(gO, k' + q'- p\ k\ /) + h.c.) , (67c) 

where the hermitian conjugation acts on the implicit flavour indices in Q. In order to 
extract the coefficients ^fl(c)'5^>'5R^ linearise in the chemical potentials for (. and R. 
Since Q is odd under a simultaneous exchange of £'i_3 ■(-)■ i?2,4 and 0, £•(-)■ A, R, we can 
do this calculation for //^ and directly infer the result for /xr. To linear order in 

jl3h^hfj.i = Qh^hqi , (68a) 



(e/3^1 - l)(e/3^2 _ i)(e/3S3 + 1) (q^e, + 1)2' 
fii =6l3^qi ■ (68b) 
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Using these relations, we can extract the reduced interaction rates: 

X sign(p°) sign(fc°) sign(g°) sign(A:° + g° - + g° - k^,p°) 

=7.72x 10-%F^(72y, 

where we have defined 7vert) ^^e universal value of the phase space integrals in (1^ . 

Now, for diagrams (A), (B) and (D), (E), it is possible to put three propagators on 
shell without cutting through the gauge boson propagator. Those cuts correspond to the 
interference between a Higgs and a gauge boson mediated process, and do not contribute 
to the equilibration of fiavours. Indeed, after linearising in deviations from equilibrium, 
these terms cancel in the collision term. 

Following the calculation of diagram (C), the relevant part of the self energy di- 
agram (A) then evaluates to (still, we suppress the superscript (0) on the tree-level 
propagators) 



■X 



(27r)4 (27r) 

iS>{k)iSj{p - q) -2^-i\ iA> {q)^^i{p - k)iA>{p - q - k) X h , (70) 

where we used that iAa^u = —Qfiui^A- To facilitate the comparison with the result 
for (C), we can use that in thermal equilibrium (with vanishing chemical potentials), 
iA^(/c) = iA'^(— fc), and iAA{k) = iA^{k). Inserting this into the collision term and 
taking the trace over Dirac indices, we find 

Cf'^' = ibl + YeY,gl) J ^^2vr5(^>^)2vr5(P)2vr5(g^) (71) 

x2n6{{k + q — p)^)sign(p°) sign(A;°) sign(g°) sign(A;° + g° — p°) 
x^(g°,A:° + g°-/,A:°,/), 
where we have used that 

tr [Pr{2P -2]^-f,{f- 0f\ ={p- qf{p - kf for {p - k - qf = ^ = p" = q^ = Q , 

(72) 

which again cancels with the denominators of the off-shell propagators. The contribution 
from diagram (B) can easily be seen to be identical to that of diagram (A). Comparing 
with Eq. fl67al) . it follows that 

^fi(A).. ^ ^fi(B).. ^ _^f,(A).R _ ^fl(B).R ^ _ (1^2 ^ ^y^y^^T^ ^flO^ _ (73) 



25 



Another way of verifying that these contributions come with the same numerical co- 
efficient is to note that diagrams (C) and (A), (B) descend from three-loop vacuum 
diagrams, that are identical in terms of spin 0, 1/2 and 1 propagators (up to exchanges 
of £ and R and of the SU(2)l and U(l)y gauge bosons). These vacuum diagrams are 
(up to the d^p integration) recovered through the integration over dp^. 

The relevant contribution from diagram (D) to the self energy is given by 

I nyi J (27r)4(27r)4 ^ ^ 

x^iSfip - k)iS>{p - k - qh.iS^ip - q)iA>{k)iA<{-q)h . 

The relevant cut through diagram (E) can be brought into the same form. When inserting 
this into the collision term, we obtain 

xtr[PR7'^(^ -^){f-]^- <l)^^{f - <l)f\ 

X p>{p - k - q)iA>{k)iA>{q)hiS<{p) - (>^<)] . 

The Dirac trace is now different than for the diagrams (A), (B) and (C), but with the 
on-shell conditions, it reduces to 



tT[PnY{p -mp-^- i)li.{f - i)f\ = 2(p - kYip - qY , (76) 
such that the collision term simplifies to 

Cfi(D,E) ^ _ 2Y,YW^h' I [iS^iP - k - q)iA>Ak)iA>{q)hiS<{p) - (>^<)] 

(77) 

= - 2Y,Ynglh^ J [^Sm^^lik + q - p)iA>{q)hiS<{p) - {>^<)] , 

where in the last term, we have shifted the momentum k ^ p — k — q, and iS^ is defined 
by iS{p) = ipS{p). Comparing with the calculations for (C) and (A,B), it is then easy 
to verify that 

^fi(D)5£ ^ ^msi ^ _ ^fi(D)5R _ ^fi(E)^R ^ _^YnY,ghm^ . (78) 
As the sum of the contributions from the various diagrams, we obtain 

.fl(A+B+C+D+E)5£ _ _ ^ fl(A+B+C+D+E)<5R _ ^fl /7q\ 
/ ~ / ~ /vertex I ' 



--gl + (2FrF^ - 21110 - 4Fl>r)^?? 



Tvert 



26 



Finally we can insert the weak hypercharges in order to obtain the flavour relaxation 
rate, and find 

iLte. = - ilgl + U) X 7.72 X 10-^ X T (80) 
= -7.72 X lO-"^ X GT - 2 X 7.72 x 10"^ x g^T . 



5 1^2 Processes 

The production rate of the singlet neutrinos also includes 1^2 processes [311 [32]. In the 
limit where Mn ^ T, the processes N and i(f)* — are the main contributions 

to the production rate. At higher temperatures relative to Mjy, which is relevant for the 
weak washout regime, the thermal masses of the lepton and the Higgs bosons, 

m', = ^{3gl + gl)T\ (81) 
lo 

ml = ^{Sgl + gl + Ah] + ^\)T^ , (82) 

are of importance. These masses are understood to be effective masses valid for modes 
of momenta larger than gT . 

Following Eqn. ( 125|) . the tree level rates can be obtained from 

tr[^S-^(p)] = -\Y f ^ ^ [X,(a;,+) -X,(a;,_))] , (83) 



where 



^i±=^\Mlj + m]-ml\ (84) 



-2Ml 



(po' - Ml) {M% + mj + ml- 2Mlm] - 2mjml - 2Mlml) 



We take the singlet neutrino to be on shell, = Mfj, and 



2i(w£) 



ilogfe/'b°-.)-lUilog(e/'-^ + l) for > + 

1 log (^e^(p°+'^«) - l) - ^ log (e^^* + 1) for m^> Mj^ + mi 

otherwise 

(85) 

However, there are also processes involving the multiple coUinear emission of soft 
gauge bosons from the lepton and the Higgs boson propagators. While kinematically 
not possible in the vacuum (at LO, when all scattering particles are massless), here they 
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occur because of the thermal masses of the gauge bosons. Note that in the scattering 
diagrams calculated in the preceding Sections, we have approximated the gauge boson 
masses by zero, which means that the coUinear processes are not readily included. 

The summation of the collinear processes is derived and discussed in Ref. ^T\. Here, 
we just quote the procedure. First, we solve the integral equations 



/(Pq± 
^^^C(q^) [/(p^,p||, A;||) - /(p^ - q^, = 2pj 



ie{k\\,p)ip{p±,p\i,kii 
where 



(2vr)2 
{2ny 



C(q±) [^{p±,p\\,kii) - ^(p^ - q±,p\\,ki 



(86a) 
(86b) 



£(/C||,pj 



2p||(P|| - kii 



k\\)M% - k\\{p\\ - k\\)mi - k\\p\\ml^ 



(87) 



and 



C(qj 



/3 



82/1 



4 



92 



qi qi 



m 



D2 



1 2M 



qi qi 



The Debye masses are mD2 = %9qT and mny = -^gfT . The solution to Eqs. (186|) 
is best performed in impact parameter space and it is not straightforward. We refer to 
Ref. [21] for the details. 

Then, the production rate for singlet neutrinos is given by 



tr 



^it^{k) 



-2Y' 



1 



|k| J (27r)3 |k| - p\\ e^Pii + 1 eP^\^\-n) - 1 
xRe 



(89) 



Ikl 

■ /(P±,P|1, |k|) + -^i){p^,p\\, |k|) 



where we take k^ > 0. The integration over d^p is performed in the limits 

Ikl 



P\\± 



(90) 



N"''i 



while the integration over |p_l| is subsequently performed from zero to infinity. 

An analytic calculation or approximation of the 1 ■(-)■ 2 rates contributing to the 
production of N is presently not available, and the numerical evaluation following above 
procedure is yet time-consuming. Moreover, a numerical generalisation to the situation 
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2.5 3.0 3.5 4.0 



Figure 7: Right-handed neutrino production rate Vj^/iY'^T'^) for gi^2 gi,2, = 0, 
varying gi^2, and all remaining couplings as given in Tabled] for the scale lO^GeV (solid 
blue line). We also indicate the fit ( 19T]) (dashed red line) and indicate the rate in the 
SM (red star). 
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Figure 8: Right handed neutrino production rate Tf^/lY'^T'^) for M^r = 0, varying gi 
and all remaining couplings as given in Table [1] for the scale 10^ GeV (solid blue line). 
The dashed red line indicates the result for gi = 0, the dotted green line the fit (l9T]l . 
The rate in the SM is indicated by a red star. 

where also the external fermion line can radiate a gauge boson has not yet been per- 
formed. In order to obtain an estimate of the uncertainty in the flavour relaxation rate 
due to the radiation from the external fermion line, but also in order to investigate the 
effect of the scale {i.e. temperature) dependence of the gauge coupling constants on the 
right-handed neutrino production rate, we calculate F^r for varying values of gi and g2- 
The numerical results are presented in Figures [7] and [HI As in this work, we are inter- 
ested in the production of light (relativistic) right-handed neutrinos and the scatterings 
of massless Standard Model leptons, gauge- and Higgs-bosons, we take for the mass 
of the right-handed neutrino Mn = 0. For the Standard Model couplings implied by a 
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Higgs boson mass of 125 GeV, we thereby reproduce the value Tj^^ ^ 5 x 10~^/(y^T^) 
found in Ref. [31]. Moreover, we find that a good fit to the behaviour apparent in 
Figure [7| is provided by 

T]^^/{Y^T^) ^ 8.8 X 10"^G . (91) 

The results presented in Figure [8] indicate that this formula becomes less accurate when 
one of the couplings is very small, i.e. we observe a deviation from the relation (1911) 
when gi <^ g2- 

Regarding the flavour equilibration rate 7^, we note that the technique developed 
in Ref. [31] applies to diagrams of the type of Figure [5t^C), where the gauge radiation 
originates from internal lines of the diagram in Figure [1] only, but not directly to the 
diagrams in Figures [SI^A,B,D,E), where the gauge radiation also attaches to an external 
line. A generalisation of the methods developed in Ref. [21] is beyond the scope of the 
present work, but as we see below, the error from neglecting gauge radiation from an 
external line is quantitatively small compared to the dominating contribution from the 
t-channel exchange of fermions and the expected NLO corrections. In order to make an 
estimate for the relaxation rate of left-handed flavour 7*^, we notice that this should be 
complementary to the relaxation rate 7^ for right-handed flavour, 

/ = ^7^ (92) 

where the factor 1/2 is due to the multiplicity of left-handed leptons and Higgs bosons. 
Now, the right handed leptons R are different from the right-handed singlet neutrinos 
in that they have the weak hypercharge —1 and that within their self-energy diagram, 
no charge-conjugated particles are running. We may therefore approximate 

tr[^Sr^(p)]^ltr[^Sr(p)], (93) 

and from Figures [7] and [HI we can estimate the relative uncertainty due to inaccurately 
neglecting the U(l) weak hypercharge interactions by 15%. As it turns out that the 1 -H- 2 
rates are small compared to the sum of the 2 o 2 rates (which are in turn dominated 
by the t-channel fermion exchange), and in view of the uncertainty from neglecting 
contributions oc G^logG*"^ at the following order (NLO), it is therefore quantitatively 
sufficient to make the estimate 

tT[fT.r\p)]^\tT[fT.'r{p)]. (94) 

Substitution into Eqs. (ITTl) and ( l2Ti) and numerical evaluation then yields (for the values 
of the couplings given in Table [1] at the scale of lO^GeV) 

7L2 = 9.9x 10-^T. (95) 

When expressed as a fit similar to relation (pTI) . this becomes 

7L2 = l-7x 10-=^GT. (96) 
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Figure 9: Individual contributions and total flavour equilibration rate /T as function 
of the temperature T. The temperature dependence emerges through the renormalisa- 
tion group equation (RGE) evolution of the coupling constants which are renormalised 
at the scale 2'kT. Earlier times appear to the right (higher temperature) in this plot. 
Contributions from gauge interactions are separated into linear and log enhanced scat- 
terings as well as 1 — )■ 2 processes. The gl and gliXoggi"^) terms are included in the 
corresponding solid and dotted blue lines and in the total rate. 

6 Phenomenological Implications 

6.1 Flavoured Leptogenesis 

The full LO flavour equilibration rate has been calculated here for the first time. Adding 
the individual contributions, it can be expressed as 

^fl = ^fl{0)^^ + ^ms^ + ^fiCR)^^ + ^fl^^^^^ (97) 

= 1.32 X 10"^ X h^^T + 3.72 x 10~^ x GT + 8.31 x 10"^ x G(log G-^)T 

+ 4.74 X 10"^ X glT + im x 10^^ x gl{\ogg:['^)T + 1.7 x 10"^GT, 

where G = |(3(7| + gf). From Section [3731 we have used the numerical fit (l60i) since it is 
better behaved for small values of p^. 

The running of the couplings let depend non-trivially on the temperature. In 
Figure |9] we show the individual contributions as well as the total flavour equilibration 
rate as function of the temperature, in the region of T = 10^ GeV to 10^^ GeV. Note that 
the tree-level 1 — > 2 rate is zero in this temperature regime, since the thermal masses for 
the Higgs and for the leptons leave no phase space for a decay process. Only once the 
collinear emission contributions are included a finite 7f^2 obtained. 

We see that the rate is largely dominated by the terms linear in G, which in turn 
receive their dominant contribution from the t-channel fermion exchange, as it is shown 
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in detail in Section I3l3l The total rate varies between 5.19 x 10~^T and 4.83 x 10~^T 
(depending on the renormalisation scale). This is close to some estimates that were 
previously used in the literature, namely 5 x 10~'^T in [291 Si]; but smaller than the 
value of 1.75 x 10~^T used in [TT]. In the light of the uncertainties due to the unknown 
NLO corrections, the improvement over the popular estimate of the flavour relaxation 
rate as 5 x lO'^T may not be dramatic. We emphasise however that in contrast to the 
latter number, the relaxation rate obtained here derives from a systematic LO calcu- 
lation. The estimates in Refs. [29| HI] go back to the work [16], where the t-channel 
divergences from fermion exchange are apparently regulated using the thermal fermion 
masses, but without including their widths, such that the numerical agreement is per- 
haps coincidental. Unfortunately, the estimation of the size of the NLO correction is not 
straightforward, such that presently, we cannot state a reliable number for the theoret- 
ical uncertainty, cf. the recently performed NLO calculation of the photon production 
rate in the quark-gluon plasma |l3]. If we were taking the discrepancy between the semi- 
analytical expression f lSSjl for the t-channel scattering rate and the numerical fit fIBU]) as 
an indication of the theoretical uncertainty, we would estimate it as 15%. 

When rescaling the results of Ref. [11] using the present LO value for the rate of 
flavour relaxation, we may conclude that the unflavoured description of Leptogenesis 
applies to masses for the decaying right-handed neutrino M^r ^ 10^^ GeV, whereas it 
may be treated as fully flavoured (no correlations between the r-lepton doublets and the 
remaining two doublets) below masses of ^ 10^^ GeV. This latter value lies below 
the sometimes estimated mass [15] of Mn ^ 10^^ GeV, where a flavoured description is 
assumed to be valid. Note that these are estimates obtained for two particular points in 
parameter space in Ref. [11], and we do not quote a quantitative error when applying 
either the unflavoured or flavoured description for a mass of the decaying right handed 
neutrino between 10^^ GeV and 10^^ GeV. A more systematic and quantitative study of 
the transition between the flavoured and unflavoured regimes will be performed in the 
future. 

The strong dependence on the gauge coupling strength also implies some model de- 
pendence of the flavour equilibration rate. In many extensions of the SM, the running of 
the gauge couplings is modified, such that the charged lepton Yukawa interactions will 
equilibrate at different temperature scales. A detailed phenomenological study of such 
variants of flavoured Leptogenesis will be presented elsewhere. 

6.2 Results for the Right-Handed Neutrino Production Rate 

The integrated rate for the production of right-handed singlet neutrinos can be obtained 
from the results of Sections [3] and H] by performing the integration with a different 
weight, as explained in Section 12.51 The main purpose of performing this integration 
is to illustrate the size of the corrections compared with the tree-level calculation, and 
to facilitate comparison with the results of Ref. [32] • For a precise numerical study of 
Leptogenesis, one should instead use the differential distribution function /Ar(|p|), since 
different momentum modes are equilibrated on different time scales, which can lead to 
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a modification of tlie resulting lepton asymmetry [T9] . 

Tfie contributions to tlie neutrino production rate from Higgs and fermion mediated 
scatterings, from vertex-type diagrams and from collinearly enhanced 1^2 processes 



are respectively given by 

= 3.15 X 10-^ X GT^ + 5.22 x 10"^ x hlT^ , (98) 

^N{i) ^ 2.24 X 10"^ X GT^ + 4.14 x 10"^ x G{\ogG-^)T^ , (99) 

7^rt = 3.15 X 10"^ X GT^ , (100) 

= 8.8 X 10"^ X GT^ , (101) 



where as before, G = + \g'l- The rate for is calculated using the analytical 

decomposition into linear and logarithmic contribution analogous to the procedure ex- 
plained in Section 13.31 If instead one extracts the coefficients from a direct integration 
of (jH]) and by performing a numerical fit, one obtains 

^N{i) = 1,57 X 10-3 X GT^ + 3.67 x 10"^ x G log G'^T^ . (102) 

See the end of Section [331 for a detailed discussion of how these two methods and results 
compare. Our numerical results should also be compared with those obtained in Ref . [32] . 
While we agree with the coefficient of the terms proportional to the top quark Yukawa 
and of the logarithmic term, we obtain a significantly larger coefficient for the term 
linear in G. Summing the terms linear in G from the contributions (p8l) - fll00p . we obtain 
2.87 X 10-3 X GT^ compared to 1.00 x lO^^ x GT^ in Ref. [32]. Using the resummed result 
for 7^(^) instead the linear coefficient is 2.20 x 10"^ x GT^, still deviating significantly 
from Ref. [32] • In summary, our result for the production rate of right-handed neutrinos 
is 

7^ = 3.08 X 10-3 X GT^ + 3.67 x 10-^ x G log G'^T^ + 5.22 x 10-^ x h^T^ . (103) 

In Figure [TOl we show the different contributions to production for vanishing 
density of A^, fN{p) = 0. We see that while the coUinear enhancement of production 
is significant when compared to the tree level rates, the gauge mediated scatterings 
dominate at high temperatures {z -C 1), while the tree level (f) — )■ Ni decay plays a 
negligible role for N production. Before z > 1, the M^v = approximation used for 
calculating the scattering rates becomes invalid. Results for the 0{g'^) corrections for N 
production in this non-relativistic regime were recently reported in Refs. [381 139] . 

Since the tree level rates are small compared to the scattering contributions in the 
high temperature regime, a precise calculation of the lepton asymmetry in this regime 
not only needs to take scatterings into account for the washout effects, but also should 
include them in the calculation of the GP asymmetry. This is particularly important 
for Leptogenesis in the weak washout regime and should therefore be addressed in the 
future. 
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Figure 10: Individual contributions and total production rate for zero N 

density, as function oi z = M^/T and Mjv = 10^^ GeV. The horizontal lines show the 
individual and total contributions from 2 — )■ 2 scattering processes. The green solid 
(dashed) lines show the tree level rate for 1—7-2 processes with (without) thermal 
masses for the Higgs and lepton doublet, while the thick green line shows the 1 — )■ 2 
rate including the collinear enhancements. Note that as before, the 2 — )■ 2 scattering 
rates are calculated in the = approximation and therefore only valid for z < 1. 
Different from Figure [H time flows from left to right here as the Universe is cooling. 

7 Conclusions 

In this work, we have presented an approach based on the OTP formalism to calculate 
interaction rates at finite temperature from 2PI self energies. This we have employed to 
calculate the flavour relaxation and the right-handed neutrino production rates relevant 
for Leptogenesis scenarios. The main results of the present paper are: 

• We have shown that finite temperature interaction rates can be calculated pertur- 
batively in the 2PI formalism. Using this approach, the t-channel divergences in 
diagrams with fermion exchange are automatically regulated. The linear and log- 
arithmic dependencies on the gauge-coupling square introduced by such processes 
can be extracted both analytically and by a numerical fit from our calculation. 

• The LO results for the flavour relaxation rates and right-handed neutrino produc- 
tion rates and their dependence on the gauge and top Yukawa couplings are given 
in Section El These expressions can easily be used for obtaining interaction rates in 
other models, since the dependencies on the temperature, gauge-couphng evolution 
and hypercharge assignment are given explicitly. 

• We find that both rates are largely dominated by the fermion-mediated t-channel 
scatterings, and are therefore sensitive to the RGE evolution of the gauge couplings. 
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The coefficients vary by about 15% depending on whether they are extracted ana- 
lytically or numerically. Since the numerical fit fl60|) includes an incomplete account 
of higher order effects through resummation, this deviation can be taken as an in- 
dication for the magnitude of higher order effects. 

• While a detailed study of the impact on flavour effects remains to be done, we 
can already conjecture in combination with the work [11] that the unfiavoured 
description of Leptogenesis will be valid for as low as 5 x 10^^ GeV, somewhat 
lower than previously estimated. 

The production rate of right-handed singlet neutrinos for vanishing density was pre- 
viously calculated in Ref. [52]. While we agree with the coefficients of the GlogG and 

terms, we find a much larger coefficient for the contributions linear in G which dom- 
inates the overall production rate. This also reduces the absolute importance of the 
coUinear scatterings that were first shown to be relevant in [31]. Further investigations 
are required to identify the origin of this discrepancy. To this end, we note that our 
semi-analytical extraction of the terms linear and logarithmic in G agrees within the 
expected accuracy with a numerical evaluation of the phase-space integrals, which can 
be taken consistency check. 

Our calculations have been performed in the approximation that all masses of the 
external particles are vanishingly small. This is a good approximation for the flavour 
relaxation rate, where all particles have masses 0{gT) that only lead to small, NLO, 
corrections since interactions in the plasma are dominated by hard processes with 0{T) 
momentum exchange. For the right-handed neutrino production rate, this approximation 
breaks down when > T, i.e. at times where z > 1. In this regime 1 — )■ 3 decays 
of are kinematically allowed and exhibit divergences when the emitted gauge boson 
is soft or coUinear with one of the other decay products. These divergences should be 
cancelled by virtual corrections that arise when only two propagators are put on shell 
in the diagrams in Figure O For the T -C M/v limit, this cancellation was recently 
demonstrated in [38| 139] . Using the approach presented here, it has been shown that the 
cancellation of soft and coUinear divergences occurs for any temperature [IQ]. 

Extensions and applications of the calculation of LO scattering and production rates 
for light (compared to the temperature) particles include a calculation of the CP- 
violating rate in the weak washout regime of Leptogenesis. Besides, a systematic and 
quantitative study of the transition regime between flavoured and unfiavoured Leptoge- 
nesis may now be performed in combination with the description of flavour decoherence 
developed in Ref. [11]. The results of this work may therefore serve as a basis for LO 
calculations of the baryon asymmetry of the Universe in scenarios where so far, only 
estimates have been available. 
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A RGE evolution of couplings 

The gauge and top quark Yukawa couplings at temperature scales relevant for Leptoge- 
nesis differ significantly from their values at the Electroweak scale. The RGE evolution 
is well known in the SM. At the one-loop level, the gauge couplings at a scale are given 
by 

«r'(/^) = «r'(Mz) - ^ log , (104) 

where hi = 41/10, 62 = —19/6 and 63 = —7, and the aj(M^) are the values of the 
couplings at the Electroweak scale. The couplings are defined as cti = ^gf/^Air) and 
Q^2,3 = fl'i 3/(471"); following the conventions used in the context of Grand Unification. 

The RGE equations for the top Yukawa and the Higgs quartic coupling at the one 
loop level are given by |47] 

(106) 

where we suppress the dependence of the couplings on the RGE scale /i. These are 
evaluated numerically using the following input values for the gauge and top Yukawa 
couphngs at the Electroweak scale g8]: ai{Mz) = 0.0169, 0^2 (M^) = 0.0338, a^{Mz) = 
0.1184 and ht{Mz) = 0.998. 

The Higgs quartic coupling depends on the Higgs boson mass through A = {rrih/v)'^ /2, 
where v = 174 GeV. Current experimental constraints indicate that = 125 GeV, 
which we take as the value of the Higgs boson mass for the rest of this analysis. 

The evolution of the couplings up to scales /i = 10^^ GeV is shown in Figure [TTJ 
The value of the Higgs quartic coupling gets negative at intermediate scales, which could 
jeopardise the stability of the Electroweak vacuum. A more detailed study, including two 
loop corrections and threshold effects, was presented in [H]. They find that while the 
Higgs coupling indeed becomes negative at high scales, it stays above the meta-stability 
bound for rrih ~ 125 GeV, and the upper limit on the reheating temperature is consistent 
with the Leptogenesis scenario. 

The values of the couplings for fi = 10^ GeV and /i = 10^^ GeV are given in Table [TJ 
These values are used for the numerical analysis in Section 16.11 The running of the cou- 
plings in the relevant temperature regime is relatively slow, such that it is not necessary 
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Figure 11: One-loop evolution of coupling constants in the SM. Note that gi = 
■\/3/5 (7i^GUT, such that gauge coupling unification at /i ~ 10^^ GeV is not evident in 
this figure. 
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Table 1: Values of the relevant SM coupling constants at scales important for Leptoge- 
nesis. 



to evolve the couplings along with the temperature in a numerical analysis. The induced 
uncertainty is of higher order in the coupling expansion. 



B Feynman Rules 

For completeness, here we list the Feynman rules that are employed to calculate the 
one and two-loop self energies directly in Wigner space. First, note that the standard 
definition, in Wigner space, is such that 




where the blob denotes the sum of all one particle irreducible (IPI) diagrams, the mo- 
mentum fiows in the direction of the arrows, and the external legs are understood to be 
amputated. 
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For fermions (and similarly for complex scalars), the two-point correlation- functions 
on the closed time path (CTP) are defined as 

iS{u,v) = {Tc{^P{u)i,{v))), (107) 

where Tc denotes time ordering along the CTP. Since i/j creates a fermion state and 
annihilates an anti-fermion, an arrow indicating particle flow will point from v to u. 
Correspondingly, we obtain the following Feynman rules for the propagators in momen- 
tum space: 

a , b =iA^\k), (108) 

« < b =iS''''{k). (109) 

Our convention that the SU(2) doublets are defined as = , (f)'^)'^ and ii = (z/j,e~)^ 
implies that the arrows point in the direction of positive (negative) electric charge flow 
for Higgs bosons (SM leptons). The momentum flows in the direction of the arrow. Since 
the gauge boson propagators are neutral, the propagators have no well defined charge 
flow. Note that we suppress all SU (2) indices in this paper. 

Interactions - Interaction vertices are derived from the Lagrangian as in conventional 
field theory. In addition, each internal vertex obtains a sign a G {+, — } that is summed 
over. The following vertices are relevant for the calculations in this paper: 



I 




where we have suppressed gauge and flavour indices. For gauge boson interactions we 
make the gauge indices explicit: 
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igiYiY^PL + YnPK)5fg for = B , 
mYPLTf, for y^ = PF^, 



- ig^Tfg {p + kf for = . 

Here are the generators of SU(2) in the fundamental representation, p, k, q denote the 
momenta and /, g are SU(2) indices. The scalar-scalar-vector (SSV) vertex is particularly 
sensitive to sign errors. The rule presented above is valid provided that both the momenta 
and the positive charge flows in the direction indicated by the arrows. 
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